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Recent  advances  in  robot  technology  have  revolutionized  the  concept  of 
manufacturing  and  assembly  systems.   These  advances  have  created  the  need 
for  new  mathematical  models  to  reflect  the  capabilities  of  the  new  techno- 
logies.  In  this  paper,  we  focus  on  the  system  design  problem  by  defining 
a  work  station  selection  and  task  assignment  problem  for  automated  assembly 
systems.   The  problem  is  formulated  as  a  zero-one  integer  program  and  a 
procedure  for  seeking  lower  and  upper  bounds  to  the  optimal  value  of  the 
Integer  program  is  described.   We  show  that  the  upper  bound  provides  a 
feasible  solution  to  the  integer  formulation  and  that  the  lov;er  bound  is 
tighter  than  the  standard  linear  programming  relaxation  of  the  integer 
formulation.   Computational  results  indicate  that  the  proposed  bounds  are 
extremely  tight.   In  fact,  in  each  of  the  42  test  problems  evaluated,  the 
lower  and  upper  bound  coincided,  indicating  that  the  optimal  solution  to 
the  integer  program  had  been  obtained. 
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With  recent  advances  in  manufacturing  and  computer  technology,  the  notion 
of  the  automated  factory  has  evolved  from  a  science  fiction  fantasy  to  a 
reality.   Some  recent  references  in  the  popular  press  which  describe  these 
advances  are  Bylinsky  (1979),  Donlan  (1980),  Friedrich,  et  al.  (1980), 
U.S.  News  and  World  Report  (1979),  Nicholson,  et  al .  (1979),  Business  Week 
(1980),  Schefter  (1980).   As  these  articles  indicate,  a  wide  variety  of  auto- 
mated equipment  for  part  fabrication,  assembly,  and  material  handling  is  now 
available  for  manufacturing  and  assembly  systems.   For  instance,  robots  are 
being  widely  used  in  welding,  spray-painting,  and  material-handling  applica- 
tions, while  other  types  of  automation  such  as  numerically-controlled  machining 
centers  are  performing  part  fabrication.   In  addition,  robots  are  expected  to 
work  along  side  assembly  workers  in  many  future  assembly  systems.   These 
technological  advances  have  created  a  new  set  of  opportunities  for  management 
scientists,  with  regard  to  both  the  design  and  the  operational  control  of 
automated  systems. 

This  paper  focuses  on  the  system  design  problem  for  automated  assembly 
systems.   An  assembly  system  performs  a  set  of  distinct  tasks  for  the  assembly 
of  a  product  family  which  typically  consists  of  a  primary  product  and  its 
model  variations.   To  perform  these  tasks,  the  system  contains  a  set  of  work 
stations  linked  together  by  a  transport  mechanism  for  moving  parts  between 
stations.   In  addition,  the  assembly  system  includes  the  detailed  specifi- 
cation of  how  a  product's  assembly  flows  from  station  to  station.   With  these 
elements  of  an  assembly  system  in  mind,  we  see  that  decisions  concerning  the 
design  of  an  assembly  system  may  be  divided  into  the  following  three  components; 
(a)   The  determination  of  the  assembly  requirements  to  be  satisfied  by  the 

system.   This  includes  specifying  the  product  family  to  be  assembled 

as  well  as  identifying  and  sequencing  the  required  assembly  tasks. 
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(b)  The  selection  of  work  stations  and  assignment  of  tasks  to  these  stations. 
Here  the  designer  not  only  determines  the  type  and  number  of  work 
stations,  but  also  assigns  the  tasks  to  these  stations. 

(c)  The  selection  of  transport  mechanisms  and  determination  of  the  physical 
layout.   The  designer  must  find  a  feasible  layout  of  the  work  stations, 
which  includes  the  choice  of  a  material  handling  system  for  transporting 
parts  between  stations  and  a  routing  plan  for  each  part  to  be  assembled. 

This  taxonomy  of  the  system  design  process  suggests  a  hierarchical  approach 
in  which  the  resolution  of  (a)  is  the  primary  input  to  (b),  while  the  deci- 
sions from  (b)  constrain  the  choices  in  (c) .   Since  the  effects  of  the  three 
decision  levels  are  interrelated,  the  decision-making  process  would  typically 
be  implemented  in  an  iterative  fashion. 

In  this  paper,  V7e  concentrate  on  (b)  ,  the  work  station  selection  and 
task  assignment  problem.   The  recent  advances  in  automated  assembly  technology 
introduce  new  complexities  and  opportunities  not  addressed  by  the  more  tradi- 
tional assembly  design  models  such  as  the  assembly  line  balancing  problem 
(Ignall  [1965]).   In  the  line  balancing  problem  we  assume  that  all  work 
stations  are  identical,  and  are  to  be  laid  out  sequentially  along  a  line.   We 
then  assign  assembly  tasks  to  these  stations  so  as  both  to  satisfy  all  prece- 
dence relationships  between  the  tasks  and  to  balance  the  total  work  load  at 
the  stations.   In  contrast,  in  the  station  selection  and  task  assignment 
problem  we  select  work  stations  from  a  set  of  nonidentical  candidate  work 
stations  and  simultaneously  assign  the  assembly  tasks  to  these  selected 
stations  to  achieve  a  prespecified  production  rate  at  lainimum  cost.   We  do 
not  presume  a  specific  layout  of  these  work  stations  so  that  considerable 
freedom  to  exchange  product  parts  between  stations  in  the  assembly  system  is 
permitted.   In  addition,  for  automated  assembly  equipment,  the  capability  for 
tool- changing  is  a  distinguishing  characteristic.   Consequently  we  explicitly 


-3- 

include  tiool  changing  in  the  problem  forinulation.   A  more  comprehensive  treat- 
ment of  our  work  may  be  found  in  Lraaar  (1980)  which  if.   an  extension  of  the 
station  selection  and  task  assignment  problem  given  in  Graves  and  Whitney 
(1979).   Future  work  will  address  parts  (a)  and  (c)  of  the  system  design 
process. 

The  remainder  of  this  paper  is  organized  into  five  sections.   The  next 
section  formulates  the  V7ork  station  selection  and  task  assignment  problem  as 
a  zero-one  integer  program.   We  illustrate  this  formulation  with  an  example 
based  on  the  assembly  of  an  automobile  alternator.    An  approximate  solution 
procedure  for  this  integer  program  is  presented  in  Sections  2  and  3.   This 
solution  procedure  consists  of  solving  tv.'o  relaxations  to  the  integer  program 
to  obtain  lower  bounds  on  the  optimal  value.   These  relaxations  are  also 
used  to  find  near-optimal  feasible  solutions  to  the  integer  program  and 
upper  bounds  on  its  optimal  value.   In  Section  4  we  report  our  computational 
experience  that  shows  that  the  approximate  procedure  is  extremely  effective 
in  finding  the  optimal  solution  to  the  integer  program.   Finally,  in  Section 
5  we  discuss  possible  model  extensions  and  refinements. 


1 .   rRESENT'ATION  OF  THE  MODEL 

.  In  this  section  we  formulate  a  zero-one  integer  linear  program  that  selects 
work  stations  and  assigns  tasks  to  these  work  stations  to  minimize  total  cost  in 
an  assembly  system  which  must  satisfy  a  preset  production  rate,  lie   then  illus- 
trate the  model  with  the  example  of  an  automobile  alternator  assembly,  and  re- 
fomiulate  the  model  as  a  generalized  programming  problem. 

Notation  and  Formulation 

Let  index  set  I  =  {1,...,M}  denote  the  M  possible  work  stations  available 
for  assembly  operations  and  let  index  set  J  =  {1,...,N}  denote  the  N  distinct 
tasks  that  must  be  performed  by  the  assembly  system.   These  N  tasks  may 
represent  a  single  product  or  subassembly,  or  possibly  a  family  of  products 
or  subassemblies.   We  define  task  0  as  a  dummy  "start-up"  task  and  task  N+1 
as  a  dummy  "completion"  task.   We  assume  that  there  is  a  strict  ordering  of 
the  tasks  such  that  task  j  must  precede  task  k  iff  j  <  k.   As  a  consequence, 
suppose  that  task  1,  2,  and  4  are  assigned  to  the  first  work  station  and 
task  3  is  assigned  to  the  second  work  station  and  must  precede  task  A.   This 
implies  that  once  the  first  work  station  performs  tasks  1  and  2,  the  product 
must  be  unloaded  and  trasferred  to  the  second  work  station  for  task  3  to  be 
performed.   Then  the  product  is  returned  to  the  first  work  station  to  complete 
task  A. 

Our  formulation  of  the  assembly  system  deisgn  problem  will  select  a 
subset  from  the  set  of  work  stations  I  and  assign  the  tasks  from  set  J  to  these 
work  stations  so  as  to  minimize  some  cost  criterion  while  achieving  a  preset 
production  rate.   This  formulation  requires  four  sets  of  parameters.   First, 
for  each  work  station  i  (id),  we  specify  an  annualized  fixed  cost  f .  represent- 
ing the  capital  investment  necessary  to  cover  the  work  station  and  its  required 
tooling.   Second,  we  specify  the  amount  of  operating  time  b. 


available  at  work  station  i  over  an  annual  period.   Typically,  b.  would 
represent  the  expected  "up-time"  for  the  work  station.   Third,  we  define 
t.  .,  (id,  0  <^  3  <  k  <^  N+1)  to  be  the  amount  of  operating  time  necessary  to 
satisfy  the  annual  requirements  for  task  j  if  work  station  i  performs  task  j 
and  task  k  is  the  next  task  performed  at  work  station  i.   The  parameter  t  ., 
typically  consists  of  three  components  -  (i)  the  actual  assembly  time  asso- 
ciated with  task  j,  (ii)-  a  load/unload  time  required  for  transferring  the 
product  to  another  work  station  if  k  does  not  equal  j+1,  and  (iii)  a  tool 
change  time  if  tasks  j  and  k  require  distinct  tools  at  work  station  i.   By 
convention,  t.„,  specifies  the  "startup"  time  at  work  station  i  when  task  k 

is  the  first  task  performed  in  each  cycle  at  that  work  station,  while  t.  .  .,,, 

1 » J  >  N+1 

gives  the  "completion"  time  when  task  j  is  the  last  task  performed  at  work 

station  i.   Finally,  similar  to  the  definition  of  t..,  ,  we  define  c... 

ijk'  ijk 

(iel,  0  <  j  <  k  £  N+1)  to  be  the  annual  operating  cost  for  performing  task  j 
at  work  station  i  when  task  k  is  the  next  task  performed  at  work  station  i. 
This  cost  includes  the  variable  assembly  cost,  the  variable  handling  cost, 
and  the  variable  tooling  cost. 

If  tasks  j,k,...,J?,  are  assigned  to  work  station  i,  the  formulation, 
for  purposes  of  calculating  the  annual  operating  time  and  cost  at  station  i, 
assumes  that  the  work  schedule  at  station  i  is  cyclic  and  follows  the  simple 
sequence  j ,k, . . . ,£, j ,k, . . . .   To  ensure  the  feasibility  of  this  operating 
sequence,  we  assume  that  the  system  carries  some  minimal  level  of  work-in- 
process  inventory.   In  practice,  it  is  possible  to  sequence 
the  work  at  each  station  so  as  to  minimize  operating  time  or  cost. 
Hence,  our  parameters  should  be  thought  of  as  conservative  approximations 
to  the  actual  costs  and  times. 

We  define  two  sets  of  decision  variables  as  follows: 
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/I  If  task  j  is  performed  at  work  r.tation  i  and  task  k  is  the  next 

)  task  to  be  performed  at  work  station  i  (id,  0  £  j  <  k  £  N+1) 

VO  othen<7ise. 

(1  if  work  station  i  is  included  in  the  assembly  design  configuration 

I  (iel) 

^0  otherwise. 


The  integer  programming  formulation  for  the  assembly  system  design  problem 

is 

M  N       N+1 

Problem  P:  Min     E      (f.y.   +     Z  Z        ex.)  (la) 

i=l        ^   "-        j=0  k=j+l      ^^^^  ^^^ 


subject   to 

M       N+1 

-         ^       X  =     1  VjeJ  (lb) 

i=l  k=j+l     ^^^ 

N       K+1 

E         E       t..,x..,      <     b.  Viel  (Ic) 

j=0  k=j+l     ^^^  ^^^^     =       ^ 

N  j-1 

E       x         -     E     x  =     0  Viel,   jeJ  (Id) 

k=j+l     ^-^'^       £=0     "-^^ 


J,   -iOk  -  ^i     =     0  ^i^I  (le) 

k=l 


N 


^     >^.-   .  M4.1   -  y.     =    0  Viel  (If) 


j=o 


iJ,N+l 


x...g{0,1}  Viel,  0<j<k<N+l     (Ig) 

y^G{0,l}  ViEl  (Ih) 

Problem  P  minimizes  the  total  annual  system  cost  given  the  annual  production 
volume  requirements.   The  objective  function  (la)  represents  the  annual  capital 
and  operating  costs  as  a  function  of  the  decision  variables.   The  first  set  of 
constraints  (lb)  specifies  that  each  task  j  be  assigned  to  a  work  station,  while 


constraints  (Ic)  ensure  that  the  amount  of  work  assij^ned  to  a  particular  station 
does  not  exceed  the  time  available  at  that  station.   Constraints  (Id)  through 
(If)  are  conservation  of  flow  constraints  that  enforce  the  precedence  relation- 
ships of  the  tasks  at  each  work  station;  thus,  if  task  j  is  assigned  to  a  work 
station  i,  then  there  exists  an  immediate  predecessor  task  I   and  successor 
task  k  to  task  j  on  station  i  (which  may  be  the  dummy  tasks  ^=0  or  k=N+l) . 
Furthermore,  these  constraints  also  relate  the  task  assignment  variables  (x..,) 
with  the  work  station  selection  variables  Cy . ) ;  if  y.  equals  zero,  so  that 
station  i  is  not  included  in  the  system  design,  then  no  task  assignments  are 
made  to*' station  i  and  the  x...  variables  are  all  zero. 

Tn  further-  hiphlipht  the  nature  of  formulation  ^,  we  note  that  by 
representing  the  work  stations  as  vehicles  and  the  tasks  as  delivery  points, 
we  may  interpret  formulation  P  as  a  special  type  of  vehicle  routing  problem 
(Dantzig  and  Ramser  [1959]).   Specifically,  nonhoraogenous  vehicles  are 
selected  to  make  deliveries  to  nodes  that  are  indexed  such  that  they  must 
be  visited  in  increasing  order  of  index.   As  described  below,  this  precedence 
structure  appears  to  make  the  problem  easier  to  solve  than  the  general 
vehicle  routing  problem  and  suggests  a  specialized  solution  procedure. 

Example 

In  order  to  illustrate  the  use  of  formulation  P,  we  apply  the  model  to 
the  assembly  of  an  automobile  alternator  (Nevins  and  \>niitney  [1978]),  v;hich 
is  depicted  in  the  exploded  view  of  Figure  1.   The  set  of  tasks  that  must 
be  performed  in  order  to  assemble  the  alternator  is  described  in  Table  1. 
The  precedence  relationships  of  the  tasks  are  prescribed  by  the  task  number; 
i.e.,  the  tasks  must  be  performed  sequentially.   Table  1  also  specifies  the 
tool  requirements  for  each  task.   A  set  of  five  work  stations,  each  being 
a  computer-controlled  robot,  is  available  to  perform  the  tasks  listed  in  the 
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Rotor  nut 


Figure  1:   Exploded  View  of  Automobile  Alternator  (Nevins  and  I-Jhitney  [1978]) 


Table  1:   Description  of  Tasks  and  Tool  Requirements  for  Automobile  Alternator 
Example 


Number  of 
Task. 

Description  of  Task 

Tool  Required 

1 

Insert 

rotor  nut 

Small  three-finger  gripper 

2 

Insert 

lock  washer 

Small  three-finger  gripper 

3 

Insert 

pulley 

Small  three-finger  gripper 

A 

Insert 

fan 

Small  three-finger  gripper 

5 

Insert 

fan  spacer 

Small  three-finger  gripper 

6 

Insert 

front  housing 

Small  three-finger  gripper 

7 

Insert 

bearing 

Small  three-finger  gripper 

8 

Insert 

retainer 

Small  three-finger  gripper 

9 

Insert 

&  tighten  3  screws 

Special  screwdriver 

10 

Insert 

rotor  spacer 

Sm.all  three-finger  gripper 

11 

Insert 

rotor 

Contracting  collet 

12 

Tighten  rotor  nut 

Rotor-nut  tightener 

13 

Insert 

rear  housing 

Large  three-finger  gripper 
Bolt  driver ^^^ 

lA 

Insert 

&  tighten  3  through  bolts 

(a)  Work  station  S.  uses  the  special  screwdriver  to  perform  task  14. 
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table.   These  robots  are  of  two  types,  general  purpose  and  special  purpose. 
The  general  purpose  robots  (referred  to  as  G   and  G^)  are  able  to  perform 
most  or  all  of  the  required  tasks.   They  have  the  ability  to  use  multiple 
tools  and  may  be  reprogrammed  to  meet  changes  in  the  production  environment. 
They  differ  both  in  their  capability  to  perform  certain  tasks  and  in  their 
costs.   In  contrast,  the  three  special  purpose  robots  (referred  to  as  S^ ,  S„, 
and  S„),  though  less  expensive  than  the  general  purpose  robots,  have  limited 
flexibility  and  are  capable  of  performing  only  a  set  of  closely  related  tasks 
using  a  single  tool.   S   is  a  special  screwdriver  station,  while  S^  and  S, 
are  pick-and-place  robots  that  are  specifically  configured  to  do  a  subset 
of  the  tasks.   See  Graves  and  VJhitney  (1979)  for  a  fuller  description  of  this 
equipment. 

We  assume  that  each  of  the  work  stations  has  120,000  minutes  of  avail- 
able time  annually,  which  are  specified  as  the  b.  parameters  in  the  model. 
If  a  work  station  is  selected,  its  time  is  consumed  by  performing  assembly 
operations,  loading  and  unloading  the  product,  and  changing  tools.   These 
three  time  components  comprise  the  t...  parameter  in  formulation  P.   Naturally. 
an  assembly  time  is  incurred  at  a  work  station  whenever  that  station  performs 
a  task.   Since  the  precedence  relationships  of  the  tasks  are  sequential,  a 
load/unload  time  is  incurred  at  a  work  station  only  when  that  station  performs 
the  first,  task  or  when  it  perfonns  a  task  j  and  does  not  perform  task  j+1. 
A  tool-change  time  is  incurred  at  a  work  station  only  when  that  station 
performs  a  task  which  requires  a  tool  different  from  the  tool  required  for 
the  succeeding  task  performed  at  that  work  station. 

Table  2  lists  the  annual  assembly,  load/unload,  and  tool-change  times 
for  the  automobile  alternator  under  consideration,  assuming  a  production  rate 
of  one  unit  per  minute  for  120,000  minutes  per  year.   We  note  that  for  robots 
G,  apd  G„  the  tool  changing  time  is  comparable  to  task  processing  time.   This 
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Table  2:  Time  Data  for  Automobile  Alternator  Example 


Time 
Component 

Work  Station  and  Annual  Time  (1000 's  of 

minutes)  ^""^ 

^1 

^2 

^1 

^2 

^3 

Assembly  Time ^ 

Task  1 

12 

8 

- 

6 

- 

Task  2 

.12 

8 

- 

6 

- 

Task  3 

12 

8 

- 

6 

- 

Task  4 

12 

8 

-  ■ 

6 

- 

Task  5 

12 

8 

- 

6 

" 

Task  6 

14 

12 

- 

" 

6 

Task  7 

12 

8 

- 

6 

Task  8 

12 

8 

-■ 

- 

6 

Task  9 

8 

8 

6 

- 

- 

Task  10 

12 

8 

- 

- 

6 

Task  11 

12 

- 

- 

- 

- 

Task  12 

12 

- 

- 

- 

- 

Task  13 

12 

- 

- 

- 

- 

Task  14 

8 

8 

6 

- 

- 

Load/Unload  Time 

2 

2 

2 

2 

2 

Tool-Change  Time^'^ 

12 

8 

- 

- 

(a)  Based  on  a  production  rate  of  one  unit  per  minute  for  120,000  minutes 
per  year. 

(b)  A  dashed  (-)  entry  indicates  an  infeasible  task/work  station  combination. 

(c)  A  dashed  (-)  entry  indicates  that  a  work  station  is  unable  to  change  tools. 
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is  a  realistic  feature  for  thesn  assembly  robots  and  illustrates  the  importance 
of  modeling  tool  changes  in  the  system  desij'n  problem.   Although  the  timing 
data  for  the  load  and  tool-change  operations  in  the  table  are  itemized  only 
by  work  station,  formulation  P  allows  this  data  to  be  different  for  each  work 
station  and  task. 

In  addition  to  timing  informations  formulation  P  requires  cost  informa- 
tion.  This  data  is  supplied  in  Table  3.   The  annual  fixed  costs  correspond 
to  the  f.  parameters  in  P  and  represent  the  procurement  of  equipment,  including 
necessary  tooling.   The  annual  variable  costs  correspond  to  the  c. .,  parameters 
in  P  and  reflect  costs  such  as  electricity  5  labor,  and  accessory  costs  v;hich 
depend  upon  the  volume  of  production.   In  Table  3  tliese  costs  again  assume  a 
production  rate  of  one  unit  per  minute  for  120,000  minutes  per  year.   Also, 
as  with  the  timing  data,  although  the  variable  costs  are  itemized  only  by 
the  work  station,  this  information  may  be  made  specific  to  both  the  work 
station  and  the  task. 

Using  the  data  supplied  above,  the  optimal  solution  to  problem  P  is 
obtained  using  the  techniques  described  in  the  following  sections  of  this 
paper.   The  optimal  total  annual  cost  for  this  problem  is  $106,300  and  the 
optimal  work  station  selections  and  task  assignments  are  displayed  in  Table  4. 
According  to  this  solution,  the  first  five  tasks  should  be  performed  on 
special  purpose  robot  S„.   The  alternator  assembly  should  then  be  moved  to 
spcial  purpose  robot  S„  vjhich  performs  tasks  6  through  8.   For  task  9,  the 
assembly  should  be  moved  to  the  general  purpose  robot  G^ .   Task  10  is 
performed  by  reloading  the  assembly  onto  robot  S„  and  then  the  assembly 
moves  back  to  robot  G   for  tasks  11  through  14  to  be  performed.   The  switch- 
ing of  the  product  between  robots  G^  and  S„  for  tasks  9,  10,  and  11  occurs, 
in  part,  because  G  does  not  have  sufficient  capacity  to  perform  tasks  9 
through  14.   Table  4  also  shows  the  percent  of  time  that  each  of  the  selected 
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Table  3:   Cost  Data,  for  Avitomobile  Alternator  Example 


Cost 
Component 

Work  Station  and  Annual  Cost  ($) 

S 

S 

^1 

^2 

s 

Fixed  Cost^^^ 
Variable  Cost^^^ 

72000 
500 

A6000 
500 

15000 
100 

15000 
200 

15000 
200 

(a)  Based  on  a  one-year  payback  period. 

(b)  Based  on  a  production  rate  of  one  unit  per  minute 
for  120,000  minutes  per  year. 


Table  4:   Optimal  Solution  for  Automobile  Alternator  Example 


Work 
Station 

Selected  in 
Optimal  Solution  ' 

Task 

Assignments  in  Optimal  Solution 

Work  Station 

Utilization 

(percent) 

^1 

yes 

9,  11,  12,  13,  14 

87 

S 

no 

- 

- 

^1 

no 

- 

- 

^2 

yes 

1,  2,  3,  4,  5 

27 

S 

yes 

6,  7,  8,  10 

23 
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work  .stations  will  be  utilized  during  an  annual  period.   Wiereas  G   is  nearly 
fully  utilized,  both  S   and  S   are  utilized  less  than  30%.   Nevertheless,  this 
is  cheaper  than  any  alternative  solution. 

We  note  that  the  above  example  is  purely  illustrative;  indeed,  it  is 
quite  simple  and  might  be  solved  by  inspection.   In  addition,  we  point  out 
that  P  does  not  explicitly  consider  the  physical  layout  of  the  work  stations. 
Rather,  in  this  model,  we  assume  that  a  feasible  layout  of  the  work  stations 
vrith  appropriate  part  feeding  devices  and  adequate  buffer  stocks  may  be 
constructed  from  the  solution  suggested  by  P.   Accordingly,  the  solution 
to  P  is  intended  as  a  useful  aid  for  the  assembly  system  designer. 

Reformula tion  of  the  Problem 

Problem  P  is  a  very  large  pure  integer  program.   Rather  than  attempting  to 
solve  P  directly,  we  instead  reformulate  this  problem  as  a  generalized  linear 
program  (Magnanti,  Shapiro,  and  Uagner  [1976])  with  binary  variables.   We  then 
solve  two  relaxations  of  the  refonnulation  -  a  pure  linear  relaxation  and  a 
mixed  integer-linear  relaxation.   In  the  succeeding  sections  of  this  paper 
we  show  that  these  relaxations  provide  not  only  very  tight  lower  bounds  to  P 
but  also  feasible  candidate  solutions  to  P  that  are  near  optimal. 

To  recast  formulation  P  into  a  generalized  programming  framework,  we 
consider  constraints  (lb)  as  the  "complicating"  constraints  and  the  remaining 
constraints  as  the  "easy"  constraints.   We  observe  that  the  "easy"  constraints 
separate  by  index  i  into  M  subprobleras.   The  feasible  region  for  subproblem  i 
is  given  by  the  set 

where  x  denotes  a  vector  with  decision  variable  x  ,  as  the  [l+j+k(k-l) /2] 

—X  IJK 

element  of  x^  fo^  ^  =  J  '^  ^  i  H+l.      Since  each  set  X  is  bounded  and  constraints 


-15- 

(Ig)  and  (l]i)  require  that  x.  and  y   be  integer,  it  follows  that  each  X^ 
contains  a  finite  set  of  points.   Hence,  by  letting  R^  denote  the  cardinality 
of  set  X.  for  iel,  we  may  express  X.  cquivalently  as 

\     =  {(2ii.yi)  :  re{l,...,R.})  Viel      (2b) 

where  (x!",y^)  is  the  r^^  point  in  set  X..   We  let  x^.,  denote  the  [l+j+k(k-l) /2: 

~-l   X  X  Xj  K 

element  of  vector  x.    for  0  ^  j  <  k  <^  N+1. 


equals  zero,  then  work  station  1  is  not  used  and  x.  is  the  zero  vector:  other- 

•    —X  ' 

wise,  for  y.  equal  to  one,  vector  x^.  may  be  interpreted  as  a  feasible  assignment 
of  tasks  to  work  station  1.   The  cost  associated  with  this  assignment  is  given  by 

N   N+1 

d.     =   J:    ^   ^iik'^iik   •  ^^^ 

"      j=0  k=j+l  ^^^  ""^^ 

We  also  define 

N+1 
a.  .   =    L       X.  ., 
•   iJ^     k=j+l  ^'^ 

as  a  zero-one  parameter  to  indicate  whether  or  not  task  j  is  included  in  the 

r   feasible  assignment  of  tasks  to  work  station  i.   With  this  notation,  the 

generalized  linear  programming  reformulation  of  problem  P  hecgpies. 

M  ^i 

Problem  Q:        Min  E   (f  y  +  I     ^ir^ir^  ^^^^ 

1=1         r=l 


subject  to 

M   ^i 

I       E  a. .  X.   =1  VjeJ  (Ab) 


i=l  r=l 


ijr  ir 


i 
^  ^^^  -  y^  =  0  Viel  (Ac) 


X^^e{0,l}  Viel,  re{l,...,R^}  (Ad) 


yj^e{0,l}  Viel  (Ae) 
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wherc  decision  variable  X.      is  a  zero-one  variable  denoting  whether  or  not  the 
ir 

r^  task  assignment  given  by  vector  >:•  is  performed  at  work  station  i. 

Formulations  P  and  Q  are  equivalent  since  there  is  a  one-to-one  correspond- 
ence between  their  feasible  solutions,  and,  since,  by  (3),  the  objective  function 
values  for  corresponding  feasible  solutions  in  the  tv.'o  formulations  are  equal. 
Moreover,  given  a  feasible  solution  to  Q,  we  may  construct  the  corresponding 
feasible  solution  to  P,  and  vice  versa. 
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2  I   PUKE  LINEAR  RELAXATION 

Like  problem  P,  problem  Q  Is  a  very  large  pure  integer  formulation. 
Accordingly,  in  this  section  and  the  succeeding  section,  we  consider  two 
relaxations  of  problem  Q.   In  this  section,  v;e  describe  the  first  relaxation, 
referred  to  as  problem  Q  ,  in  which  we  replace  the  integrality  contraints 
(Ad)  and  (Ae) ,  respectively,  with 

0  <  X^^  <  1  Viel,  red R^}      (5a) 

0  <  y.  <   1  Viel  (5b) 

Our  interest  in  problem  Q  is  threefold.   First,  by  design,  this  problem  is  a 
pure  linear  program  whose  solution  procedure  fits  naturally  into  a  generalized 
programming  algorithm  (i.e.  Dantzig-Wolf e  decomposition).   Second,  the  optimal 
value  of  Q  provides  a  tighter  lower  bound  to  the  optimal  value  of  P  than 
does  the  linear  programming  relaxation  of  P.   Third,  Q  is  used  to  seek  a 
feasible  solution  to  the  pure  integer  problem  P,  and  hence  an  upper  bound 
to  the  optimal  value  of  P. 


Solution  Procedure 

To  solve  Q  ,  we  consider  a  column  generation  procedure  (Dantzig  and 
Wolfe  [1961]),  in  V7hich  we  require  an  efficient  means  of  evaluating  each  of 
the  M  subproblems.   For  each  subproblem  i,  the  feasible  region  is  given  by 
the  set  X.  defined  in  (2);  the  objective  function  is 
N   N+1 

where  -u  =    (tt.)  is  the  vector  of  dual  multipliers  associated  with  the  constraint 
set  (Ab),   and  ir^  =   0.   Thus,  subproblem  i  for  iel  is 
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Problcra  S.:  Mln  v.(Yr) 


subject   to  (x.,y   )eX. 


We  observe  that,  for  y.  =  0  we  have  x.  =  0^  and  v.(tt)  =  0.   When  y.  =  1,  S. 
is  an  acyclic  shortest  path  problem  together  with  a  single  time  availability 
constraint.   One  approach  to  solving  problems  of  this  structure  has  been 
suggested  by  Handler  and  Zang  (1978),  who  maximize  a  Lagrangean  relaxation  of 
S.  in  which  the  complicating  constraint  (Ic)  is  relaxed;  they  resolve  duality 
gaps  by  branch- and-bound.   We  consider  an  alternative  approach  which,  for  each 
index  i,  solves  the  subproblem  with  constraints  (Ic)  through  (Ih)  as  an 
expanded  acyclic  shortest  path  probleBi  and  hence  guarantees  an  integer  solution. 

To  solve  S.  for  y.  =1,  we  construct  a  network  G,  with  node  set  N.  and 
arc  set  A..   Node  set  N.  contains  ordered  pairs  of  elements  (j,^)  in  vjhich 
the  first  component  represents  a  task  assigned  to  station  i  and  the  second 
component  is  the  cumu3  ative  time  consumed  at  work  station  i.   Symbolically, 

Nj  =  {(j,i^)  :  jc{0,...,N+l},  £e{0,...,b.}}.  (7) 

The  arcs  contained  in  the  arc  set  A.  are  of  the  generic  form  [(j,?-)>  (k,m)], 
corresponding  to  the  performance  of  task  j  at  v7ork  station  i  at  time  i   followed 
by  the  performance  of  task  k  at  work  station  i  at  time  m.   The  strict  precedence 
relationships  of  the  tasks  are  maintained  in  network' G.  by  requiring  that 

(  >  j   if  je{0,...,N} 
k  (8a) 

(  =  j   if  j=N+l 

The  time  consumed  at  work  station  i  by  performing  task  j  when  task  k  is  the 
next  task  performed  at  that  work  station   is  represented  in  network  G  by 
letting 
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[l  +   t,.,   if  je{0,..,,N} 
.   m  =       ^J^  (8b) 

\l  +  1      if  j  =  N+1. 

Hence  the  arc  [  (j  ,Ji) ,  (k,iii)  ]e:A.  if  and  only  if  (j  ,«-)  ,  (k,m)cN^  and  (j,£-),(k,m) 

satisfy  (8).   A  unit  flow  on  arc  [  (j  ,?0  ,  (k,in)  ]  corresponds  to  setting 

X,.,  =  1;  thus,  the  arc-cost  for  arc  [  (j  .i!-)  ,  (5<,in)  ]  is  (c    -  it  )  . 
ijk  iji^    J 

Using  this  structure,  each  path  from  node  (0,0)  to  node  (N+l,b^)  in 
network  G.  corresponds  to  a  solution  to  S .  in  which  y.  =1.   Stated  another 
V7ay,  each  path  represents  an  integral  assignment  of  tasks  to  work  station  i 
which  satisfies  both  the  precedence  relationships  (Id)  through  (If)  and  the 
time  availability  constraints  (Ic)  at  work  station  i.   Hence,  we  solve  S^  with 
y  =  1  by  finding  the  shortest  path  on  G..   The  shortest  path  problem  may  be 
readily  solved  by  simultaneously  constructing  and  evaluating  network  G  using 

a  forward  dynamic  programming  procedure.   Assuming  that  the  parameters  b 

2 
and  t...  are  integer,  this  procedure  has  a  computational  complexity  of  0(N  'b^) 

operations  (Lawler  [1976]).   If  the  optimal  value  of  this  shortest  path  problem 

is  less  than  -f . ,  then  these  values  are  the  optimal  solution  to  S.;  otherwise  the 

optimal  solution  to  S.  is  y.  =0  and  x.  =  0^. 

Lower  Bound 

We  wish  to  show  that  problem  0  provides  a  tighter  lower  bound  to  the 
optimal  value  of  the  integer  problem  P  than  the  linear  programming  relaxation 
of  P.   For  notational  convenience,  we  let  P  denote  the  linear  programming 
relaxation  of  P  and  let  v(A)  denote  the  optimal  value  of  any  problem  A.   Using 
this  notation,  we  claim  that 

v(P)   <  v(Q-'-)   <  v(P)  (9) 

The  right  inequality  of  (9)  is  trivial  to  obtain  since,  as  we  have  observed. 
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probleins  P  aiid  Q  are  equivalent  and  problem  Q  is  a  relaxation  of  Q. 

To  demonstrate  the  validity  of  the  left  inequality  of  (9),  we  consider 
the  Lagrangean  relaxation  (Shapiro  [1979])  L(ji_)  of  P  by  dualizing  constraints 
(lb)  to  obtain 

M  N   N+1  N 

Problem  L(TT):     Min  I      [f-Y.  +  ^    ^    ^^iik~''i^''iilJ  "^  ^  ^    ^^^^ 
1=1   ^  ^   j=0  k=j+l  ^-"^     ^     ^^  j=l  ' 

subject  to       (Ic)  through  (Ih), 

vhcre  tt  ='  (u.)  is  the  vector  of  Lagrangean  multipliers  of  constraint  (lb) 
with  Ti  ''H  0.   The  optimal  value  of  ]l   is  obtained  by  solving,  the  dual  nroblem 

Problem  D:       Max  v[L(t[)]. 
H 

As  shown  in  Magnanti,  Shapiro,  and  Wagner  (1976),  we  may  re\-nrite  problem  D 
as  a  linear  program  whose  dual  problem  is  given  by  Q  .   Thus,  v;e  have 
v(D)  =  v(Q  ).   In  addition,  we  note  that  the  linear  programming  relaxation  of 
L(2i_)  is,  in  general,  not  equal  to  L(ji_).   This  implies  that  Geoffrion's  inte- 
grality property  (1974)  does  not  hold  so  that  v(P)  ^  v(D) .   Therefore,  by 
simple  substitution,  we  obtain  v(P)  <^  v(Q  ).   Moreover,  computational  results 
show  that  v(P)  <  v(Q  )  is  possible. 


Upper  Bound 

The  gen.-^ralized  programming  procedure  used  to  solve  Q  may  also  be  employed 

to  seek  an  upper  bound  to  the  optimal  value  of  problem  P.   Specifically,  in 

each  iteration  of  the  generalized  programming  procedure,  we  tost  whether  or  not 

the  decision  variables  {X.  }  are  integer.   If  so,  then  the  solution  is  a  feasi- 
ir  ' 

ble  solution  to  P  and  the  corresponding  objective  function  value  is  an  upper 
bound  to  v(P) . 

If,  in  the  process  of  solving  Q  ,  we  do  not  find  an  integer  solution  at 
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any  iteration,  we  may  employ  a  heuristic  procedure  to  obtain  a  solution  feasi- 
ble in  P.   Assume  that  the  optimal  solution  to  Q   contains  a  fractional  task 
assignment;  that  is  0  <  X .   <  1  for  at  least  one  i  and  r.   To  resolve  the 
fractional  task  assignments  and  construct  a  feasible  solution  to  P,  we  find 
the  largest  fractional  variable,  say  Xt",  and  round  it  up  to  one;  we  reset  Xj 
to  zero  for  all  r  -/  r.   This  results  in  making  the  bask  assignment  correspond- 
ing to  vector  x-,    to  work  station  i.   We  then  select  the  next  largest  fractional 
variable  and  repeat  the  process  until  either  all  tasks  are  assigned  to  some 
work  station  or  all  available  work  stations  have  been  selected.   In  the  former 
case,  i'f  a  task  is  assigned  to  more  than  one  work  station,  then  that  task  is 
reassigned  solely  to  the  work  station  mth  minimal  operating  cost;  in  this 
manner  a  feasible  solution  is  easily  found.   In  the  latter  case,  this  heuristic 
has  failed  to  find  a  feasible  solution;  but  this  failure  is  not  a  serious 
worry  since  we  can  always  find  a  feasible  solution  by  expanding  the  candidate 
set  I  to  permit  more  work  stations. 

If  a  feasible  solution  to  P  is  obtained,  then  its  objective  function 
value  may  be  compared  to  the  lower  bound  for  v(P)  to  obtain  a  conservative 
a  posteriori  measure  of  the  near-optimality  of  the  candidate  solution.   Compu- 
tational experience  has  indicated  that  not  only  do  these  upper  bound  procedures 
nearly  alv;ays  generate  a  feasible  solution  to  P,  but  also  that  the  objective 
function  values  of  these  solutions  are  very  close  to  the  optimal  value  of  P. 
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3.   MIXED  n'rrF.GER-I.TNEAIl  RE].A.\'7\TI0N 

In  the  preceding  section  we  demonstrated  that  problem  Q  ,  the  linear 
programming  relaxation  of  problem  Q,  may  be  used  to  seek  lower  and  upper  bounds 

to  v(P),  the  optimal  value  of  the  integer  program  P.   In  this  section,  we 

2 

consider  a  second  relaxation,  referred  to  as  problem  Q  ,  in  which  we  linearize 

the  A.  _  variables  but  maintain  the  integrality  requirements  for  the  y.  variables 

2 
used  to  select  v7ork  stations.   In  other  words,  problem  Q  is  problem  Q  with 

constraints  (4d)  replaced  with  constraints  (5a) .   This  second  formulation 

merits  investigation  since  it  provides  a  tighter  lower  bound  to  v(P)  than  does 

Q  ,  and-.its  solution  permits  us  to  generate  additional  candidate  solutions  to  P. 

2 

Of  course,  since  Q"  is  a  mixed  integer-linear  formulation,  its  solution  requires 

a  greater  computational  effort  than  that  for  0  .   Below,  vje  outline  the  solution 

2 
procedure  employed  to  solve  0   and  describe  the  method  of  determining  lower 

and  upper  bounds  to  v(P).   The  procedure  for  seeking  the  upper  bound  also 

supplies  a  feasible  candidate  solution  to  problem  P. 


Solution  Procedure 

2 
Problem  Q  has  M  integer  variables,  namely  the  y.  variables  for  iel.   Thus, 

we  may  solve  this  relaxation  by  employing  a  standard  branch-and-bound  procedure 

(Garfinkel  and  Nemhauser  [1972])  on  these  variables.   The  binary  enumeration  tree 

for  the  branch-and-bound  procedure  consists  of  M  levels  in  which  each  level  specifies 

a  y,  variable  as  either  zero  or  one.   For  notational  convenience,  we  define  the  sets 
\ 
Iq  =  (i  :  y^  -  0,  iel}  , 

Using  this  notation,  a  node  in  the  enumeration  tree  may  be  specified  by  the 

pair  (1^,1  )  where  I„,  I  are  disjoint  sets.   For  each  such  node,  we  define 

2 
problem  Q  \oth  constraints  (4e)  replaced  with 
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Vielp 

(11a) 

Viel^ 

(lib) 

Viel  -  (IqDT-j) 

(lie) 

'i  =  « 


y^  e  {0,1} 

-2  2 

We  let  Q  (1^,1  )  denote  the  linear  programming  relaxation  of  Q  (I  , I  )  in  which 

constraints  (lie)  are  linearized. 

The  branch-and-bound  procedure  enumerates,  either  explicitly  or  implicitly, 

all  nodes  (I  , I  )  in  which  I^UT,  =  I.   We  initiate  the  procedure  by  solving  Q 

_2 
which  is  equivalent  to  solving  Q  (I_,I  )  with  I^^  =  I^  =  0.   We  then  evaluate 

'■-2 
problem  Q  (I  , I  )  at  each  node  (I„,I^)  either  by  solving  the  linear  program 

via  the  generalized  programming  procedure  described  in  Section  3  or  by  solving 

its  dual  problem  via  a  subgradient  optimization  of  a  Lagrangean  relaxation 

(Shapiro  [1979]).  In  the  latter  case,  we  evaluate  the  Lagrangean  relaxation  given  . 

by  (10)  in  which  the  integer  restriction  on  the  y^  variables  in  (Ih)  are  replaced 

with  (11).   Based  on  the  value  of  Q^d^.I^),  the  branch-and-bound  procedure  may 

_2 
fathom  a  node  (and  all  of  its  descendants)  in  the  enumeration  tree  if  vlQ  (Iq^I^)] 

2   . 

exceeds  the  value  of  the  best  knovm  candidate  solution  to  Q  ;  in  addition, 

-2  2 

a  node  is  also  fathomed  if  the  solution  to  Q  (1,1)  is  feasible  in  Q  ,  which 

yields  a  new  candidate  solution.   If  the  procedure  cannot  fathom  a  node  (1^,1^^) 

by  solving  Q^(I  ,1  ),  then  the  procedure  "branches"  on  a  y^  variable  in  which 

i(5l„UI,  to  form  two  new  nodes  -  one  in  which  y.  is  one  and  the  other  in  which 

_2 
y  Is  zero.   We  note  that  when  IqUI,  =  X,  the  solution  to  Q  (Iq.I-j^)  always 

-2  2 

fathoms  node  (Iq.Ij^)  since  Q  (Iq,I]^)  =  Q  (^q'^I^  ^"^  ^°  ^"^  solution  to 

-9  '       2 

Q  (I  ,1  )  is  also  feasible  in  Q  .   The  branch-and-bound  procedure  continues 

until  all  nodes  are  fathomed,  at  which  point  the  best  candidate  solution 

2 

obtained  is  the  optimal  solution  to  Q  . 

_2 
In  solving  Q  (I„,ip,the  choice  between  using  generalized  linear 
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progrnmniing  versus  solving  a  dual  problem  via  a  subgradient  optimization  proce- 
dure  depends  on  the  desire  for  generating  feasible  candidate  solutions  to  the 
integer  program  P.   The  subgradient  optimization  procedure  is  very  efficient 
and  provides  a  lov.^er  bound  to  Q  (I  ,1^)  at  each  iteration;  hov;ever,  since  it 
solves  the  dual  problem,  this  procedure  yields  a  feasible  solution  to 
Q  (I„,I  )  only  when  the  optimal  solution  is  found.   On  the  other  hand,  the 

generalized  linear  programming  procedure  is  a  primal  procedure,  which  is 

-2 
less  efficient  but  wliich  finds  a  feasible  solution  to  Q  (I^jl^)  at  each  itera- 
tion.  These  intermediate  solutions  are  also  feasible  in  P  if  all  the  X. 

ir 

variables  are  integer.   In  addition,  the  heuristic  procedure  described  in 
Section  2  may  also  be  used  to  seek  a  feasible  solution  to  P.   Accordingly,  we 
use  the  primal  procedure  to  evaluate  Q  (I  ,1^)  wlien  I^  =  I   =  0  or  when 


method.   VJe  have  found  such  a  procedure  to  be  very  effective  in  finding  near- 
optimal  candidate  solutions  to  P. 

In  suinmary,  our  solution  procedure  for  P  is  a  heuristic  procedure  based 

12  1 

on  the  two  relaxations,  Q   and  Q  .   We  first  solve  Q  by  a  column  generation 

procedure  to  obtain  a  lower  bound  on  v(P).   In  this  process  we  generate  feasi- 
ble candidate  solutions  to  P,  the  best  of  which  provides  an  upper  bound  on 

v(P).   If  this  upper  bound  is  sufficiently  close  to  the  lower  bound,  we  stop 

2 
the  procedure.   Otherwise,  the  second  relaxation  Q  is  solved  by  a  branch-and- 

bound  procedure  to  obtain  a  better  lower  bound  on  v(P).   Again,  in  this  process, 

we  seek  feasible  solutions  to  P,  the  best  of  which  we  take  as  our  solution. 

The  near-optimality  of  the  best  feasible  solution  can  be  determined  by  compari- 

2 
.goa  with  the  lower  bound  given  by  v{Q  ).       . 
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4.   COMPIITATT  ONAL  PI' I^FOKiMANCE 

We  now  examine  the  coinpiitational  effectiveness  of  the  solution  procedure. 
To  do  this,  we  solved  a  series  of  test  problems  which  used  the  alternator 
example,  presented  in  Section  1,  as  a  prototype.   Below,  we  present  the  experi- 
mental plan  undertaken  to  generate  the  test  problems  and  report  the  computa- 
tional effort  required  in  their  solution. 

Experimental  Plan 

We  consider  a  production  setting  in  which  three  general  purpose  types  of 
work  station  and  five  special  purpose  types  of  work  station  are  available  to 
assemble  a  product  requiring  the  completion  of  up  to  20  tasks.   The  Appendix 
contains  the  cost  data,  time  data,  and  tool  requirements  for  the  test  problems. 
This  data  is  based  on  product  information  from  the  automobile  alternator 
example,  and  on  available  cost  and  performance  information  for  robots  currently 
available  to  industry.   Three  different  characteristics  -  the  number  of  multiple 
general  purpose  robots,  the  number  of  tasks,  and  the  production  rate  -  are 

varied  in  the  experimental  plan  in  order  to  evaluate  the  effectiveness  of 

1      2 
relaxations  Q  and  Q  under  alternative  problem  specifications.   Changes  in 

the  number  of  work  stations  and  tasks  vary  the  complexity  of  the  problems 

while  changes  in  the  production  rate  simulate  alternative  production  environments. 

The  first  characteristic  is  the  number  of  general  purpose  work  stations. 

It  is  often  necessary  to  allow  the  selection  of  more  than  one  work  station 

of  a  particular  type.   One  way  to  represent  multiple  work  stations  is  to 

define  a  separate  binary  y,  variable  for  each  of  the  identical  work  stations. 

Hence  to  permit  the  choice  of  up  to  seven  identical  work  stations  requires 

seven  y.  variables.   However,  an  alternative  method  for  representing  identical 

work  stations  is  to  define  multiple  work  stations  consisting  of  2   identical 

work  stations,  n=0,l,2,....   For  instance,  a  double  work  station  (n=l)  repre- 
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sents  tv,'0  Identical  work  stations,  and  has  twice  the  fixed  cost  and  twice  the 
capacity  of  a  single  station.   Hence  to  permit  the  choice  of  up  to  seven  identi- 
cal work  stations  requires  the  specification  of  only  three  multiple  work  stations; 
a  single  (n=0) ,  a  double  (n=l) ,  and  a  quadruple  (n=2)  work  station.   Thus,  the 
second  method  permits  representation  of  multiple  work  stations  with  fewer  y. 
variables. 

In  the  experimental  plan,  the  number  of  general  purpose  work  stations  is 
divided  into  three  cases.   The  first  case  allovjs  at  most  one  of  each  general 
prupose  work  stations,  so  that  M  equals  8.  'The  second  case  allows  at  most 
three  o.f  each  general  purpose  work  stations.   This  is  done  by  defining  a  double 
station  for  each  general  purpose  work  station  and  thus  M  equals  11.   The  third 
case  allov;s  at  most  seven  of  each  general  purpose  work  station  and  is  modelled 
in  a  similar  manner  so  that  M  equals  14.   The  second  characteristic  in  the 
experimental  plan,  the  number  of  tasks,  is  divided  into  five  categories  by 
considering  selected  subsets  of  the  20  tasks.   These  five  categories  contain 
two  categories  in  which  N  is  10  (tasks  1-10  and  tasks  11-20) ,  two  categories 
in  vjhich  N  is  15  (tasks  1-15  and  tasks  6-20),  and  one  category  with  N  equal 
to  20  (tasks  1-20).   The  third  characteristic,  the  production  rate,  is  set 
at  three  levels  -  2.0,  4.0,  and  6.0  units  per  minute. 

In  summary,  implementation  of  this  experimental  plan  results  in  the 
investigation  of  a  total  of  45  test  problems  ranging  in  size  from  8  to  14  work 
stations  and  from  10  to  20  tasks. 


For  each  of  the  45  test  problems,  we  solve  the  two  relaxations,  Q  and 
Q  ,  in  order  to  examine  both  the  tightness  of  the  lower  and  upper  bounds  to 
the  optimal  value  of  the  integer  formulation  P  as  well  as  the  computational 
effort  required  to  obtain  these  bounds.   Of  the  45  test  problems,  42  are 
feasible  in  formulation  P.   Table  5  reports  the  deviation  of  the  lower  and 
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Table  5:   Tightness  of  Lower  and  Upper  Bounds  Cor  Test  Problems 


Relaxation 

Absolute  Deviation  of  Bounds  fron  Optimal  Value  of  Problem  P  (%) 

Lower  Bound 

Upper  Bound 

Maximum 

Minimum 

Average 

Median 

Maximum 

Minimum 

Average 

Median 

Q^ 

20.0 

0.0 

4.1 

0.0 

15.5 

0.0 

1.5 

0.0 

Q^ 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

(a)   Based  on  42  feasible  test  problems. 


Table  6:   Computational  Performance  for  Test  Problems 


Relaxation 

Computational  Performance  Measures 

Total  Running  Time  ,  , 
(CPU  seconds)   ^^ 

Number  of  Problems  Evaluated,  . 
in  Enumeration  Tree 

Maximum 

Minimum 

Average 

Median 

Maximum 

Minimum 

Average 

Median 

Q^ 

393.0 

10.0 

59.7 

38.0  - 
40.0 

- 

~ 

- 

- 

Q^ 

1086.0 

10.0 

151.6 

68.0  - 
70.0 

198.0 

0.0 

28.5 

0.0 

(a)  Based  on  42  feasible  test  problems. 

(b)  Using  a  PRIME  400  minicomputer. 

_2 

(c)  Excludes  the  initial  linear  program  Q  (0,0) 
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1      2 

upper  bounds  obtained  by  solving  Q   and  Q   from  v(P) ,  the  optimal  value  of  P. 

For  problem  Q  ,  although  the  lower  bound  ranges  as  much  as  20  percent  below 
v(P)  and  the  upper  bound  ranges  as  much  as  16  percent  above  v(P),  these  bounds 
are  equal  in  over  one-half  of  the  feasible  test  problems.   Even  more  interesting 
is  the  fact  that  in  33  of  the  42  problems  solved  by  Q  ,  the  candidate  solution 

associated  with  the  upper  bound  is  the  optimal  solution  to  P  (as  verified  by 

2 
solving  Q  ).   Thus,  in  three-fourths  of  the  problems,  the  optimal  solution  to 

the  pure  integer  problem  P  is  obtained  simply  by  solving  the  pure  linear 

problem  Q  . 

For.  the  test  problems  in  which  the  lower  and  upper  bounds  obtained  by 

1  2 

solving  Q   do  not  coincide,  application  of  Q"  yields  a  further  refinement  of 

these  bounds.   Remarkably,  after  completing  the  branch-and-bound  procedure  used 

2 
to  solve  Q  ,  the  lower  and  upper  bounds  are  equal  for  each  of  the  kl   problems, 

and  hence  the  optimal  solution  to  P  is  obtained.   The  generalized  programming 

procedure  is  by  far  the  most  common  method  of  obtaining  the  upper  bound.   It 

produces  the  best  upper  bound  in  four-fifths  of  the  feasible  test  problems 

while  the  heuristic  procedure  produces  the  best  upper  bound  in  only  one-fifth 

of  these  problems. 

2 
Although  solving  Q  yielded  the  optimal  solution  to  P  for  all  test  problems, 

2 
examples  exist  where  solving  Q  does  not  solve  P.   In  addition  to  this  compu- 
tational experiment,  we  have  used  this  solution  procedure  for  a  few  assembly 
system  design  projects  for  actual  industrial  products.   For  these  projects  we 

have  occasionally  (roughly  once  in  every  15-20  problem  instances)  found  tliat 

2 
solving  Q  did  not  solve  P.   However,  in  all  instances  the  solution  procedure 

2 
generated  a  feasible  solution  to  P  with  a  solution  value  within  2%  of  v(Q  ), 

which  is  a  lower  bound  to  v(P). 

The  computational  effort  required  to  solve  each  of  the  A2  feasible  test 

1      2 
problems  by  Q  and  Q   is  supplied  in  Table  6.   The  first  measure  of  the  compu- 
tational effort  provided  in  the  table  is  the  running  time  in  cpu  seconds  on  a 
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TRIME  AOO  time-shared  ininicon.puter.  Althoush  cpu  f.ime  is  machine  dependent 
(1  cpu  second  on  a  PRIME  400  is  roughly  equivalent  to  0.2  cpu  seconds  on  an 
IBM  370/158)  and  nay  fluctuate  with  the  number  of  users  sharing  the  machine, 

it  still  provides  a  rough  basis  of  comparison  for  computational  results.   For 

2 
instance,  we  observe  that,  on  average,  the  running  time  of  Q   is  roughly  twice 

that  of  Q""".   Since  the  Q"  running  time  includes  the  solution  of  the  initial 

problem  Q^  (0,fl),  which  is  equivalent  to  Q  ,  we  see  that  the  computational  effort 

required  to  solve  the  initial  linear  program  is  approximately  the  same  as 

_2 
that  for  the  evaluation  of  the  remaining  Q  Ci^Q^i)    problems  in  the  enumeration  tree, 

-2 

The  second  measure  in  Table  6  reflects  the  actual  number  of  Q  (I_,I  ) 

_2 
problems  [other  than  the  initial  problem  Q  (0,0)]  that  were  solved  in  the 

branch-and-bound  procedure  used  to  evaluate  the  enumeration  tree.   Even  though 

this  number  is  as  large  as  198  for  one  test  problem  (with  14  available  work 

stations  and  20  tasks),  we  observe  that,  on  average,  the  enumeration  tree  was 

evaluated  fairly  effectively.   We  attribute  this  fact  to  the  tight  upper  bounds 

_2 
produced  by  the  initial  problem  Q  (0,0). 
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5.   MODEL  REFINEMENT 

Tlie  reported  work  on  the  station  selection  and  task  assignment  problem  in 
the  context  of  assembly  system  design  is  part  of  an  ongoing  research  effort 
concerning  the  development  of  adaptable/programmable  assembly  systems;  this 
vork  is  being  done  in  conjunction  with  a  research  team  at  the  Charles  Stark 
Draper  Laboratory  in  Cambridge,  MA.   Our  future  work  on  the  assembly  system 
design  problem  will  include  the  refinement  of  the  current  model,  given  by  P, 
to  capture  better  the  essence  of  the  v/ork  station  selection  and  task 
assigpjnent  problem.   In  addition,  we  intend  to  explore  other  aspects  of 
the  system  design  process,  as  mentioned  in  Section  1,  and  to  try  to  under- 
stand better  the  remarkable  performance  of  our  approximate  solution  procedure, 
as  reported  in  the  yjrevious  section. 

The  desire  to  refine  the  current  formulatioii  is  a  consequence  of  our 
experience  using  this  model,  and  an  earlier  version  of  it  (Graves  and  VJhitney 
[1979]),  on  industrial  studies  in  which  the  research  team  attempted  to  design 
an  assembly  system  for  an  actual  product  or  family  of  products.   I'Jhitney,  et  al, 
(1979)  and  (1981)  detail  some  of  these  case  studies.   These  studies 
indicate  that  the  current  model  is  generally  very  useful  to  the  design  engineer. 
In  particular,  the  model  helps  the  assembly  system  designer  to  consider  both 
the  economic  tradeoffs  and  the  feasibility  restrictions  inherent  in  the  design 
of  the  system.   Furthermore,  the  existence  of  the  model  as  a  design  and  evalua- 
tion tool  has  made  the  designer  more  conscious  of  the  economic  criterion, 
and  presumably  has  resulted  in  better  system  designs.   Despite  this  success, 
the  case  studies  have  also  identified  three  restrictions  on  the  current  model 
which  limit  its  effectiveness  as  a  design  tool. 

First,  the  current  formulation  does  not  accurately  model  a  work  station 
for  which  an  upper  limit  exists  on  the  number  of  different  tasks  which  can  be 
assigned  to  it.   This  restriction  is  common  for  pick-and-place  (fixed-stop) 
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robots  which  can  access  only  a  finite  number  of  locations;  in  addition  this 

restriction  can  occur  due  to  limitations  on  space  for  fixtures  and  part 

feeders,  as  well  as  limitations  on  computer  memory.   To  represent  this 

restriction  we  need  to  augment  P  with  the  following  set  of  constraints 

N   N+1 

I  I       X,  .,   <  n,  Viel  (12) 

j=0  k=j+l  ^'^     =   ^ 

where  n.  is  the  maximum -number  of  tasks  possible  at  work  station  i.   We  can 
reformulate  this  augmented  problem  as  a  generalized  program  wherein  sets 
X.  defined  in  (2)  also  incorporate  (12).   This  redefinition  of  X^  results 
in  a  more  complicated  subproblem  for  the  generalized  linear  programming 
procedure;  otherwise,  the  solution  procedure  given  In  Sections  3  and  4 
remains  unchanged. 

Second,  the  current  formulation  does  not  model  economies  of  scale,  in 
both  operation  time  and  cost,  that  occur  when  similar  tasks  are  clustered 
at  a  single  work  station.   Formulation  P  assumes  that  distinct  tasks  are 
to  be  considered  independently,  except  for  possible  tool  change  and  load 
times  associated  with  the  preceding  and  succeeding  tasks  performed  at  that 
work  station.   However,  this  need  not  be  the  case  since  a  particular  work 
station  may  be  capable  of  performing  a  cluster  of  tasks  more  efficiently 
than  is  suggested  by  the  individual  task  times.   To  model  these  economies 
of  scale  we  need  to  assume  that  the  task  clusters  are  not  too  numerous  and 
that  these  clusters  can  be  identified  a  priori.   Then  the  dynamic  programming 
procedure  uh;ed  to  solve  the  subproblems  S-  ™ay  be  modified  to  evaluate  each 
task  cluster  separately. 

Finally,  the  current  formulation  does  not  model  the  fixed  tooling  costs 
associated  with  the  purchase  and  maintenance  of  tools  necessary  for  the  system 
design.   This  cost  depends  upon  the  task  assignments  which  dictate  the  tools 
needed  at  each  work  station.   To  model  the  fixed  tooling  cost  in  P,  we  need  to 
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introduce  a  nev;  set  of  zero-oiKi  variables  to  denote  Xvihether  or  not  a  parti- 
cular tool  is  needed  at  a  v;ork  station.   We  cannot  solve  the  resulting  formula- 
tion with  the  present  solution  method.   However,  we  note  that  the  inclusion  of 
fixed'  tooling  costs  in  the  model  is  critical  only  when  these  costs  are  signi- 
ficant and  when  candidate  system  designs  may  require  that  certain  tooling  be 
duplicated  at  various  work  stations. 
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APPENDTX 

Tables  Al,  A2,  and  A3  contain,  respectively,  the  cost  data,  the  time 
data,  and  the  tool  requirements  used  in  the  experimental  plan  outlined  in 
Section  5. 


Table  Al:   Cost  Data  for  Test  Problems 


Cost 
Component 

Work  Station  and  Annual  Cost  ($) 

"^ 

^^2 

S 

\ 

=2 

h 

\ 

S 

Fixed  Cost^^^ 
Variable  Cost^^^ 

72000 
500 

A6000 
500 

70000 
500 

15000 
100 

15000 
200 

15000 
200 

15000 
200 

15000 
200 

(a)  Based  on  a  one-year  payback  period. 

(b)  Based  on  a  production  rate  of  one  unit  per  minute  for  120,000  minutes  per 
year. 
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Table  A2:   Time  Data  for  Test  ProblemE 


Time 
Component 

Work 

Statioi 

and  An 

nual  Time  (1000 's  of 

minutes 

)(a) 

^1 

^2 

s 

^1 

^2 

^3 

^4 

h 

Assembly  Time^ 

Task  1 

12 

8 

10 

- 

6 

- 

- 

- 

Task  2 

12 

8 

10 

- 

6 

- 

- 

- 

Task  3 

12 

8 

10 

6 

-- 

- 

- 

Task  4 

12 

8 

10 

- 

6 

- 

- 

- 

Task  5 

12 

8 

10 

-- 

6 

- 

- 

- 

Task  6 

14 

12 

12 

- 

- 

6 

- 

- 

Task  7 

12 

8 

10 

- 

- 

6 

- 

- 

Task  8 

12 

8 

10 

- 

- 

6 

- 

- 

Task  9 

8 

8 

10 

6 

- 

- 

6 

- 

Task  10 

12 

8 

10 

- 

- 

6 

- 

- 

Task  11 

12 

- 

10 

_ 

~ 

- 

- 

- 

Task  12 

12 

_ 

10 

- 

- 

- 

- 

6 

Task  13 

12 

- 

10 

- 

- 

- 

- 

- 

Task  14 

8 

8 

10 

6 

- 

- 

- 

- 

Task  15 

12 

8 

10 

- 

- 

_ 

6 

- 

Task  16 

12 

8 

10 

- 

- 

- 

6 

- 

Task  17 

12 

8 

10 

- 

- 

- 

6 

- 

Task  18 

12 

8 

10 

- 

- 

- 

- 

6 

Task  .19 

12 

8 

10 

- 

- 

- 

- 

6 

Task  20 

12 

8 

10 

- 

- 

- 

- 

6 

Load /Unload  Time 

2 

2 

2 

2 

2 

2 

2 

2 

Tool-Change  Time 

12 

8 

10 

- 

- 

- 

- 

- 

(a)  Based  on  a  production  rate  of  one  unit  per  minute  for  120,000  minutes 
per  year. 

(b)  A  dashed  (-)  entry  indicates  an  infeasible  task/work  station  combination. 

(c)  A  dashed  (-)  entry  indicates  that  a  work  station  is  unable  to  change  tools. 
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Table  A3:   Tool  Requirements  for  Test  Problems 


Number  of  Task 

Tool  Required 

1 

Small  three-finger  gripper 

2 

Small  three-finger  gripper 

3 

Small  three-finger  gripper 

4 

Small  three-finger  gripper 

5 

Small  tliree-finger  gripper 

6 

Small  three-finger  gripper 

7 

Small  three-finger  gripper 

8 

Small  three-finger  gripper 

9 

Special  screwdriver 

10 

Small  three-finger  gripper 

11 

Contracting  collet 

12 

Rotor-nut  tightener 

13 

Large  three-finger  gripper 

14 

Bolt  driver ^^^ 

15 

Special  screwdriver 

16 

Special  screwdriver 

17 

Special  screwdriver 

18 

Rotor-nut  tightener 

19 

Rotor-nut  tightener 

20 

Rotor-nut  tightener 

(a)  Work  Station  S^  uses  the  special  screwdriver  to  perfoi 
task  14. 
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